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ABSTRACT 

We discuss the dust distribution within photoionized regions. Assuming a geometry with a central 
dust cavity, which is strongly suggested by the literature, we can estimate the cavity's radius from the 
ratio of the infrared and radio fluxes by using a simple transfer model of Lyman continuum photons. 
We apply the method to a sample of the Galactic H II regions. The estimated typical radius of the dust 
cavity of the Galactic compact H n regions is about 30% of the Stromgren radius. Taking account of 
uncertainties both of the observational data and the model, we can reject a dust distribution lacking 
a central cavity. Therefore, the dust cavity model is supported independently of the previous works. 
We discuss the formation mechanism of such a dust cavity and its detectability by present and future 
infrared facilities. 

Subject headings: dust, extinction — H II regions — infrared: ISM — radiative transfer — 
submillimeter 



1. INTRODUCTION 

Dust grains exist everywhere. The inside of H n regions 
is not an exception (e.g., Ishida & Kawajiri 1968; Harper 
& Low 1971; Melnick 1979). The radiation from celestial 
objects is always absorbed and scattered by the grains. 
Without correction for the dust extinction, we inevitably 
underestimate the intrinsic intensity of the radiation, and 
our understanding of the physics of these objects may be 
misled. In order to evaluate the amount of the dust ex- 
tinction in detail, it is important to know the distribution 
of dust as well as its optical properties. 

Theoretically some attempts to reveal the dust distri- 
bution in H II regions have been made to date (Mathews 
1967, 1969; Gail & Sedlmayr 1979a,b). According to Gail 
& Sedlmayr (1979b), the radiation pressure on the dust 
grains by the central source causes a central hole of dust 
and gas. The dust distribution has also been observation- 
ally investigated. From the observed flux ratio of H (3 to 
continuum, O'Dell & Hubbard (1965) have shown that the 
gas-dust ratio decreases with radius in the Orion nebula. 
From radial photometric profiles at U-band of some H n 
regions, Nakano et al. (1983) also suggest that a dust de- 
pletion zone is in the central region (See also Yoshida et 
al. 1986, 1987). Moreover, the central dust cavity is sup- 
ported by fitting to the observed infrared (IR) spectral 
energy distribution (SED) of H n regions (Chini, Kriigel, 
& Kreysa 1986; Chini, Kriigel, & Wargau 1987; Church- 
well, Wolfire, & Wood 1990; Faison et al. 1998; Ghosh et 
al. 2000). Recently, such dust geometry is preferred not 
only by IR SED fitting but also by the radial photometric 
profiles at submillimeter range (Hatchell et al. 2000). In 
addition, we have some evidence of central gaseous cavities 
from radio observations (Terzian 1965; Wood & Church- 
well 1989). Since dust and gas are well coupled each other 
(Gail & Sedlmayr 1979b), the presence of a gaseous cavity 
indicates a dust cavity. 

Although the dust cavity is likely to exist, it has not 
been observed directly in the far-infrared (FIR) band, 



which is dominated by the thermal emission from dust. 
This is because the FIR observations are not advanced 
relative to other wavelengths, and moreover, the spatial re- 
solving power of FIR imaging is generally weak. Recently, 
more detailed observational data have been obtained by 
ISO (Infrared Space Observatory) and SCUBA (Submil- 
limeter Common- User Bolometer Array), and many suc- 
cessive facilities are planned in the near future (SIRTF, 
SOFIA, ASTRO-F, ALMA, etc.). By using these new 
powerful facilities, we may observe the dust distribution 
directly. We need to estimate the detectability of the cav- 
ity by these facilities. 

On the other hand, we have dveloped a method for de- 
termining the fraction of Lyman continuum (LC) photons 
contributing to hydrogen ionization (Inoue, Hirashita, & 
Kamaya 2001; Inoue 2001, hereafter Paper I and Paper II, 
respectively). In paper I and II, we show that only half of 
LC photons from the central source in an H n region ion- 
izes neutral hydrogens, and the rest are absorbed by dust 
grains within the ionized region. If a dust cavity exists 
in the ionized region, however, the efficiency of the dust 
absorption for the LC photons becomes small. In this pa- 
per, we examine whether the estimated photon fraction is 
reproduced by a model nebula with a central dust cavity. 

In the next section, we describe our method for examin- 
ing the dust distribution, especially the central dust cav- 
ity's radius, by using the fraction of LC photons contribut- 
ing to hydrogen ionization. In section 3, we apply it to 
the sample H n regions in the Galaxy. Then, we discuss 
the formation mechanism of the dust cavity and its de- 
tectability by some IR facilities in section 4. Finally, we 
summarize our conclusions in the last section. 

2. METHOD 

2.1. Photon Transfer 

We assume an ideal spherical nebula composed of only 
hydrogen and dust because the effect of helium is not so 
important for the determination of the dust distribution. 
This will be discussed in more detail in the last paragraph 
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of section 2.3.2. Moreover, we assume the volume number 
densities of hydrogen and dust to be constant, nu and n<j, 
respectively. If densities of neutral hydrogen atoms, pro- 
tons, and electrons, are represented by n H o, n p , and n e , 
respectively, we have figo = (1 — x)nu and n p = n c = xnn 
in terms of the ionization degree, x, which varies with the 
radius, r, from the central exciting source. In addition, we 
assume the frequencies of all LC photons to be that of the 
Lyman limit, that is, all quantities are treated monochro- 
matically in our transfer model. This is because we obtain 
only the number of LC photons from observations of re- 
combination lines or thermal radio emissions from H II 
regions. 

We introduce a dimensionless radius, y = r/R§, where 
Rs is the classical Stromgren radius which is defined as 

1/3 



Rs = 



(1) 



Here TVlc is the intrinsic production rate of LC photons 
of the central source, and as is the Case B recombination 
coefficient of hydrogen. 

We consider the change of the number of LC photons 
passing through the shell whose inner and outer radii are 
y and y + dy, respectively. If the number of LC photons 
entering the shell per unit time denotes N^dy), its change 
in the shell is 

dN L c(y) = -A r Lc(j/){^H°SH + (1 - Lu)n d Sd}Rsdy 

n p n e aidy , (2) 

where sh and Sd are cross sections for LC photons of neu- 
tral hydrogen atoms and dust grains, respectively, u> is the 
dust albedo for these photons, and oei is the recombination 
coefficient to the ground state of hydrogen. We neglect the 
asymmetric effect of the dust scattering for simplicity. 
Since the shell is in ionization equilibrium, 

NLc{y)niios H R s dy = AnRly 2 n v n c a A dy , (3) 

where a a is the coefficient of hydrogen recombination to 
all levels. We note that a a — olb = ct\. By equation (3) 
and replacing n H o with (1 — x)nn, equation (2) is reduced 
to 

dN LC (y) = -N LC (y){(l - x)— + XjnnSnRsdy , (4) 
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where 



(5) 



This is a correction term for existence of dust grains in the 
nebula. If we define the total (gas + dust) optical depth 
for LC photons as 

(6) 



dr = {(1 -x) h X}ny.snRsdy , 



Then equation (3) is rc- 



we obtain N^dy) — A^c e ~ 
duced to 

1 — x 4TrR 2 i y 2 nHOtA 
x 2 sh^Vlc e ~ T 

Now we define the ionization parameter as 



U = 



1 / a B 2 n l iN L c\ 
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where c is the light speed. Then equations (6) and (7) are 
reduced to 

dT 3cS » {{ l- x )^L+X}U, (9) 



dy 



and 



a B 



l-x 



aAy 2 e T 



(10) 



respectively. 

Finally, we determine two recombination coefficients of 
hydrogen, a a and as- The coefficients are a function of 
the electron temperature, T e (Spitzer 1978). Fitting values 
in Table 5.2 of Spitzer (1978), we can approximate them 
to the following equation; 



= 2.06 x 10" 



a, -logT c /4000K 



where i = 



(11) 

A or B, and 



for 2000 K < T c < 16000 K 
a a — 2 .40 and = 1.64. The uncertainties of values de- 
termined by these equation are less than 1% for the above 
range of T c . 

2.2. Dust Properties 

Let us estimate the dust effect that appears in the term, 
X. We adopt the dust extinction law in the Galactic inter- 
stellar medium (ISM) developed by Weingartner & Draine 
(2001), who have derived a size-distribution function of 
dust grains composed of carbonaceous and silicate popu- 
lations by Draine & Lee (1984) and Li & Draine (2001). 
Their carbonaceous grain includes PAH (polycyclic aro- 
matic hydrocarbon) molecules. Weingartner & Draine 
(2001) have developed the extinction laws for various val- 
ues of which is the ratio of the visual extinction to 
color excess. 

We approximate the optical depths for LC photons 
by hydrogen and dust to those at 912 A. According to 
Weingartner & Draine (2001), the dust cross sections per 

unit hydrogen nucleus at 912 A i.e. ^dSd 12A /^H are 
2.515 x 10~ 21 cm 2 , 1.676 x 10~ 21 cm 2 , and 1.087 x 10~ 21 
! for Ry = 3.1, 4.0, and 5.5, respectively. 1 We note 



cm 



-22 



that these values are calibrated by Ai/Nr = 2.6 x 10 
cm 2 , where A^h is the hydrogen column density. Also, 
the dust albedos at 912 A are 0.2451, 0.2694, and 0.2989 
for Ry = 3.1, 4.0, and 5.5, respectively. If we adopt 
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6.3 x 10 18 cm 2 , equation (5) is reduced to 
f 3.0 x 10~ 4 for R v = 3.1 



X 



(12) 



1.9 x 10~ 4 for R v = 4.0 • 
I 1.2 x 10~ 4 for R v = 5.5 
Here, we introduce a parameter for the dust distribu- 
tion, the dimensionless dust cavity's radius, j/d, because 
the dust distribution may have a radial dependence as 
strongly suggested in the literature (§1). For simplicity, 
we consider an extreme case that no dust is inside the 
cavity's radius, j/d, and the gas-dust ratio is constant out- 
side of it. That is, the dust term, X, is zero within j/d, 
and constant value shown above outside of it. 

Now, we can solve the photon transfer problem. If the 
ionization parameter, U, is given, and r = and x = 1 at 



We adopt the extinction laws with the parameter sets recommended by Weingartner & Draine (2001). For Rv = 3.1, we adopt the law with 
the parameter set of Case A and 6c = 6.0 X 10~ 5 , which represents C abundance of very small grain population. Also, we adopt the laws 
with the parameter sets of Case B and 6c = 4.0 X 10~ 5 , or Case B and 6c = 3.0 X 10~ 5 for Rv = 4.0, or 5.5, respectively. Values of other 
parameters are tabulated in Table 1 in Weingartner & Draine (2001). 
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y = are adopted as the boundary condition, we obtain 
the ionization degree, x, at any radius, y, from equations 
(9) and (10) for a set of Ry and T c . The dust inner ra- 
dius, yd, is only one free parameter in our photon transfer 
model. 

2.3. Ionizing Photon Fraction 

As investigated in Paper I and II, not all of the intrin- 
sic LC photons emitted by the central source of H n re- 
gions is consumed in hydrogen ionization. In other words, 
the number of LC photons contributing to the ionization, 
i.e. ionizing photons, is always less than the intrinsic one. 
Then, we define the ionizing photon fraction, /, as 



/ 



Kc 

A LC 



(13) 



where N^ c and TVlc are the production rates of the ioniz- 
ing and intrinsic LC photons, respectively. For the ionizing 
photons, the production rate, iV^c ^ s determined from the 
observational data, whereas the intrinsic production rate, 
iVix, cannot be obtained observationally. Let us deter- 
mine the fraction, /, via two ways below. 

2.3.1. From Transfer Model 

First, we determine the fraction from our transfer model. 
According to Paper I, the fraction, /, is reduced to yf for 
a perfectly ionized nebula (i.e., x = 1 for y < and x = 
for y > J/;), where y\ is the dimensionless radius of the 
ionized region. 

Total number of proton in a nebula is J xuuAtt R§y 2 dy in 
our transfer model. On the other hand, the proton number 
is equal to (4Tr/3)R^yfnn in the perfectly ionized nebula. 
Since our solutions of the photon transfer model indicate 
that the nebula is regarded as the perfectly ionized nebula, 
we can equate above two quantities each other. Thus, we 
obtain 



/ 



Jo 



xy dy, 



(14) 



where y max is the position to stop calculating. We can 
choose the value of y max arbitrarily if it is larger than 
In this paper, y max = 1.2 is adopted as a constant for sim- 
plicity. Therefore, we can determine / as a function of 
the dust inner radius, yd, for three Ry values, if an ion- 
ization parameter, U, is given. We note that the fraction 
determined from above equation (14) becomes equivalent 
to that in section 2 of Paper I if yd = and u) = 0. 

In Figure 1, we show such results for various sets of U 
and Ry when we adopt T c = 6600 K which is a mean elec- 
tron temperature of the sample H II regions described in 
section 3. We find that / decreases if yd decreases. This is 
because the dust optical depth is large as the cavity's ra- 
dius is small. We also find that / decreases if U increases, 
which is equivalent that / decreases if Alc or nu increases 
(eq.[8]). This indicates that the ionizing photon fraction 
for a dust-gas ratio 2 becomes smaller as the number of LC 
photons from the central source becomes larger or as the 
density of the surrounding medium becomes higher. 



The relation between / and yd is also a function of Ry, 
as shown especially in Figure 2 (a). That is, the ionizing 
photon fraction for a fixed dust cavity's radius becomes 
smaller, as the value of Ry becomes smaller. This is be- 
cause the dust extinction law of a smaller Ry shows a 
steeper rise at the far-ultraviolet range (Weingartner & 
Draine 2001). Since a higher value of Ry is suitable for a 
higher density region (Cardelli, Clayton, & Mathis 1989), 
we select Ry = 5.5 for our sample H II regions discussed 
in section 3. Moreover, we see the variation of the f-yd 
relation for the change of the adopted electron tempera- 
ture, T c . We find in Figure 2 (b) that / becomes smaller 
as T e becomes higher, because recombination coefficients 
becomes larger as T c increases. 

2.3.2. From Observational data 

Next, we determine the ionizing photon fraction, /, from 
the observational data. According to Paper I and II, we 
can estimate / from the ratio of the IR luminosity, Ltr, 
to the ionizing photon production rate, N^ c , (equation [2] 
in Paper II) unless a significant amount of LC photons es- 
capes from H II regions. 3 However, in the derivation in 
Paper II, the stellar SED are assumed to be the Planck 
function. Thus, we overestimate the intrinsic production 
rate of LC photons, and then, we underestimate / in Paper 
II. 

In this paper, we adopt a realistic stellar spectra, the 
Kurucz's ATLAS 9 spectra with solar abundance and tur- 
bulence speed of 2 km s _1 . Then, if a Salpcter's IMF 
(0.1-100 Mq) is assumed, we derive the following formula; 
0.375 + 0.625e 
~ 0.250 + 4.42Ltr 6 /A[ C49 ' ( } 

where e denotes the average efficiency of dust absorption 
for nonionizing UV (A > 912 A) photons, and -Lir,6 and 
A^ c 49 are the total IR luminosity and ionizing photon 
production rate normalized by 10 6 Lq and 10 49 s _1 , re- 
spectively. The value of / obtained from the above equa- 
tion is typically 1.2 times larger than that from the equa- 
tion (2) in Paper II because of the modification of adopted 
stellar SED. When we estimate e from the observed color 
excess (see Paper II for detail), we determine / of individ- 
ual H II region from its Lir and N^ c . It is worthwhile to 
nate that since the above equation is based on a simple 
energy conservation, the obtained / is independent of the 
dust distribution, the structure of H II regions and so on. 

Once / is obtained, we can determine the intrinsic LC 
photon production rate by using equation (13), and can 
determine the ionization parameter, U. Using the ob- 
tained U, we determine yd so as to reproduce / determined 
from equation (15) by using our photon transfer model (i.e. 
eq.[14]). In the next section, we determine a typical yd for 
a sample of the Galactic H n regions. 

We discuss the effect of the presence of helium. Of 
course, helium atoms absorb LC photons as well as hydro- 
gen and dust. However, the fraction of stellar LC photons 
ionizing helium is small, < 10%. Even if all helium ion- 
izing photons are consumed by only helium, the fraction 



of photons ionizing hydrogen, /, becomes at most 10% 

2 In this paper, we adopt a dust-gas ratio corresponding to Aj/Nn = 2.6 X 10 — 22 cm 2 (see section 2.2). Since the dust-gas ratios are different 
from one region to others (e.g., Stanimirovic et al. 2000 for the SMC), we should determine the ratio for the individual H II region. However, 
we use the above dust-gas ratio for all sample regions, because we cannot determine it individually. 

3 For details of the derivation, see also Inouc, Hirashita, & Kamaya (2000). 
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smaller than that determined by equation (15). Thus, the 
determined yd decreases at most 10%. Since not all helium 
ionizing photons are absorbed by helium and most of he- 
lium recombination photons can ionize hydrogen (Mathis 
1971), the effect of helium is negligible for our analysis. 

3. RESULTS 

In order to determine the inner radius of the dust dis- 
tribution, yd of an H II region, we need various data; the 
electron density and temperature, the ionizing photon pro- 
duction rate, the infrared luminosity, and Ry- Here, the 
value of Ry is assumed to be constant, 5.5. The effect of 
the change of Ry is discussed quantitatively in the last 
part of this section. 

We find a suitable data set of the Galactic H n regions in 
Simpson et al. (1995), which examined the physical prop- 
erties of 23 H II regions in the Galaxy observed by the 
Kuiper Airborne Observatory. We select 13 H n regions 
from the sample of Simpson et al. (1995). The adopted 
selection criterion is that the measured electron number 
density exceeds 1000 cm~ 3 . The selected regions are clas- 
sified into the 'compact' or 'ultra-compact' H n regions 
from their densities (Habing & Israel 1979). The spherical 
geometry of our transfer model may be not so bad for such 
sample, although even the 'ultra-compact' H n regions 
sometimes show non-spherical shapes (Wood & Church- 
well 1989). For more evolving H n regions than the 'com- 
pact' class, which often show very complex shapes, e.g., 
'blister', 'bubble', 'champagne', etc., the spherical assump- 
tion is not valid. The observed and measured quantities 
of the sample regions are tabulated in Table 1. 

First, we estimate the ionizing photon fraction, /, 
for each sample region. The ionizing photon pro- 
duction rate, N^ Cl is estimated from the flux den- 
sity of free-free radio emission at 5 GHz by the 
formula in Condon (1992): Aj^/s" 1 = 8.88 x 
10 46 (rc/10 4 K)-°- 45 ( l //5GHz) 01 (^/kpc) 2 (^/Jy). The 
total IR (8-1000 /im) luminosity, Ltr, is estimated from 
the observed IR (18-160 /im) flux by an upward factor of 
1.26, assuming the dust SED to be the modified Planck 
function of 30 K with the emissivity index of 1. The con- 
version factor varies between 1.03-1.26 if we change the 
adopted parameter sets of the dust temperature (30-40 
K) and the emissivity index (1-2). Thus, the conversion 
causes at most 20 % uncertainty of the derived IR lumi- 
nosity. The determined N^ c and Lir are given in columns 
(2) and (3) of Table 2, respectively. 

Since we do not have the color excess, Eb-v, of sample 
regions, we assume e = 1, which correspond to a typical 
Eb-v 01 the Galactic H n regions observed by Caplan 
et al. (2000) (See Paper II). The determined / may be 
overestimated if the true value of e is smaller than unity. 
Indeed, if we set e = 0.8, the value of / becomes about 90 
% of the determined here. Now, we can determine / for 
each sample region via equation (15), and give the values 
in column (4) of Table 2. The mean value and standard 
deviation of / for 13 sample regions are 0.55 and 0.22, re- 
spectively. Thus, only half of LC photons contributes to 
hydrogen ionization in the H n regions of the Galaxy. This 
is consistent with the results in Aannestad (1978), Paper 
I and II. 



Since the ionizing photon fraction, /, is indirect mea- 
surement value, many uncertainties of adopted quanti- 
ties accumulate. Suppose a quantity, 77, is represented by 
a functional form, <j>(x\, X2, •••), where X{ is more basic 
quantities. If uncertainties of xi are represented by Si, 
and all quantities of independent each other, the 

uncertainty of r\ can be expressed by ^2(d(j)/dxi) 2 5i. In 
this way, the uncertainty of Lir/Al C is estimated to be 
about 30 % typically. The main causes of the uncertainty 
are that of the adopted radio flux densities, about 20 % 
(Simpson et al. 1995). Taking account of the uncertainties 
of the coefficients in equation (15) and the parameter, e, 
we estimate a typical uncertainty of / to be about 33 %, 
which corresponds to ±0.2 for the mean /. Thus, a some- 
what large value of / for G298. 22-0.34, in stead of / < 1 
by the definition, may be caused by this uncertainty. 

Next, we determine the dust cavity's radius, yd, of sam- 
ple regions. By using the estimated /, we can convert the 
ionizing photon production rate, A LC , into the intrinsic 
LC one, Alc, that is, Alc = (V/)-^lc- Then, we calcu- 
late the ionization parameter, U, of the individual sample 
region by equation (8) shown in column (5) in Table 2. 
Here, we determine yd so that / estimated from equation 
(15) is equal to / calculated from our transfer model (i.e., 
eq. [14]). 4 The obtained values of yd are given in the last 
column in Table 2. The mean value of the determined yd 
and its standard deviation are 0.39 and 0.23, respectively. 
The uncertainty of the individual value of yd caused by 
that of the estimated / is about ±0.2. 

In Figure 3, we show histograms of the derived quan- 
tities. From the panels (c) and (d), we can divide the 
sample into two groups: one is the 'ultra-compact' H 11 
regions, whose Stromgren radius, (or the ionized radius, 
r;) is smaller than 0.5 pc (0.4 pc). The other is the 'com- 
pact' H 11 regions with more larger radii. Here, we should 
note that G298.22-0.34 is removed from the following dis- 
cussions because of its anomalous value of /, and SI 56 
is also removed because of its very small luminosity rel- 
ative to others. The two classifications are indicated in 
Figure 3 and Table 1. In Table 3, we summarize the mean 
properties of the sample of the two groups. Although our 
classification is based on the radius of sample regions, the 
classification is caused by the difference of the electron 
density, and are consistent with that of Habing & Israel 
(1979). 

For the 'compact' sample, we estimate the mean value 
of yd to be 0.30 ± 0.12, where the error is the sum of 
two uncertainties of the mean value; one results from the 
variation among sample regions, and the other from the 
uncertainty of individual value. The estimated error is 
coincident with the sample standard deviation. For the 
'ultra-compact' sample, the mean value of yd is estimated 
to be 0.48 ± 0.26. The error is the same mean as that of 
the 'compact' sample. 

The non-zero mean value of yd obtained above may re- 
sults from a somewhat large uncertainty of individual ob- 
ject. If we assume that the sample regions compose a nor- 
mal population, we can test whether the mean yd of the 
population is regarded to be zero. If the variance of yd for 
the population is estimated to be square of the uncertainty 
of individual value, (0.2) 2 , we can reject the hypothesis of 



4 Determining y^, wc adopt the recombination coefficients calculated from the individual electron temperature by equation (ff). 
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the mean yd = with the significance level of 0.0011% for 
the 'compact' sample, and 0.0016% for the 'ultra-compact' 
sample. Therefore, even if we take account of somewhat 
large uncertainties of the observations and the model, the 
ionizing photon fraction cannot be reproduced by a nebula 
filled with dust uniformly (i.e. yd = 0). The dust distribu- 
tion with the central cavity is supported by our analysis 
independently of the IR SED fitting or the photometric 
profile analysis. 

The cavity's radius in units of the Stromgren radius (i.e. 
yd) of the 'ultra-compact' sample is larger than that of the 
'compact' sample. The real scale of the cavity's radii of 
both samples are about 0.3 pc and 0.2 pc for 'compact' 
and 'ultra-compact', respectively. We may see an evolu- 
tionary sequence of the dust cavity from 'ultra-compact' 
to 'compact'. However, these values are agree with each 
other within their error bars, and the 'ultra-compact' sam- 
ple consists of only three regions. Thus, it is uncertain 
whether the difference between the cavity's radii of two 
subsamples is real or not. 

Finally, we demonstrate the effect of the change of Ry 
on the obtained dust cavity's radius, yd- As shown in Fig- 
ure 2 (a), the value of yd increases for a fixed ionizing 
photon fraction, /, if Ry decreases. This is shown quanti- 
tatively for three sample set in Table 4. For the 'compact' 
sample, the value of yd of Ry = 3.1 is a factor of 1.7 times 
larger than that of Ry = 5.5. 

4. DISCUSSIONS 

As described above, assuming spherical symmetry and 
uniform density distribution, we obtain a typical radius 
of the central dust cavity in the 'compact' H II regions, 
0.3 times Stromgren radius. This corresponds about 40 
% of the ionized radius and about 0.3 pc. We discuss 
the formation mechanism and the detectability of the cav- 
ity. However, we should keep in mind that since real H 
II regions, even compact or ultra-compact regions, show 
a complex structure, the radius of a real cavity may be 
different from 0.3 pc obtained by our model. 

4.1. Formation Mechanism 

Three formation mechanisms of the central dust cavity 
are naturally expected: (1) radiation pressure, (2) stellar 
wind by the central source, and (3) dust sublimation. Here 
we discuss whether the above mechanisms can produce the 
cavity and which is more effective. 

The effect of radiation pressure of the central source on 
the distribution of gas and dust in H II regions is discussed 
in detail by Gail & Sedlmayr (1979b). According to them, 
the radiation force acting on the dust grain produces the 
central cavity of dust and gas (See their Fig. 7). The radius 
of the cavity is about 20 % of the ionized radius. Although 
we cannot directly compare their cavity's radius with that 
in this paper because of the difference between the adopted 
parameter sets, the radiation pressure can produce a large 
cavity. 

Also, strong stellar wind of the massive stars must con- 
tribute to produce the central cavity. Indeed, the central 
low density cavity of gas is produced by the stellar wind 
(e.g., Comeron 1997). Since dust and gas are well coupled 
with each other (Gail & Sedlmayr 1979b), the gas cavity 

5 See http://sirtf.caltech.edu/SSC/ 



indicates the dust cavity. Unfortunately, we cannot esti- 
mate quantitatively the contribution of the stellar wind to 
form the dust cavity here. It is beyond the scope of this 
paper. We will resolve which of radiation pressure and 
stellar wind is dominant mechanism to form the cavity in 
our future work. 

How about the effect of dust sublimation? Indeed, dust 
grains may sublime in H II regions because of the strong 
radiation from the central source. However, the radius of 
the dust cavity caused by the dust sublimation is about 
10~ 4 pc (Mookerjea & Ghosh 1999). This is too small to 
produce the large cavity expected here. Thus, the dust 
sublimation is not dominant mechanism to form the dust 
cavity. 

In addition, we discuss the effect of the clumpy distri- 
bution of dust. Under the clumpy distribution, which is 
more realistic than the smooth one adopted here, we may 
reproduce the ionizing photon fraction determined from 
observational data even if the dust clumps exist in the 
central area of H n regions. However we consider that 
there is also the central cavity under the clumpy distri- 
bution. This is because the radiation pressure and stellar 
wind will blows out the dust clumps from the central area. 
We note that the expected radius of the cavity in clumpy 
medium may become smaller than the determined one in 
the previous section. Such effect of the dumpiness should 
be clarified in our future work. 

4.2. Detectability of Dust Cavity 

As stated by Wright (1973), the best observational test 
of the dust cavity model is FIR (or submillimeter) observa- 
tions with high angular resolution. Assuming a spherical 
and uniform structure, we expect that the dust cavity's 
radius is about 0.3 pc for our sample of the Galactic com- 
pact H II regions. It corresponds to about 60"/(Z)/kpc), 
where D is the distance to the region from us. Here, we 
examine the detectability of such dust cavity of 'compact' 
or 'ultra-compact' H n regions by the present and future 
IR facilities. 

The future IR observational satellite of Japan, ASTRO- 
F (IRIS; Infrared Imaging Surveyor) is planned to launch 
in 2004 by the Institute of Space and Astronautical Sci- 
ence of Japan. It will offer modest angular resolutions of 
30"-50" at 50-200 /im (Murakami 1998). We may detect 
the dust cavity of H II regions located within only about 
1 kpc from us. However, ASTRO-F will survey the whole 
of the sky, fortunately. Thus, we may detect many such 
closest cavities. 

SIRTF, which is the major upcoming IR satellite 
planned to launch in July 2002 by NASA, offers much 
higher angular resolutions of 2.5"-16" at 24-160 /im. 5 If 
we observe the emission from dust being thermally equi- 
librium with ambient gas, the super resolution mode at 
70 /im, which provides us with 4.9" angular resolution, is 
the most suitable. By this observing mode, we can detect 
the dust cavity of H II regions within about 12 kpc from 
us, i.e. almost all H II regions in our Galaxy. Therefore, 
SIRTF will resolve the dust cavity, and clarify the dust 
distribution in H II regions in detail. 

We may also detect the dust cavity at submillimeter. 
For example, SCUBA provides us with much higher angu- 
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lar resolution, 8" at 450 /xm, or 14.5" at 850 fim. Hatchell 
et al. (2000) observed the Galactic 'ultra-compact' H II 
regions by SCUBA, and present high resolution images 
of them. Hunter et al. (2000) also present high resolution 
images of the 'ultra-compact' H n regions obtained by Cal- 
tech Submillimeter Observatory. However, we cannot find 
the evidence of the cavity from their images. This may 
be because the sample regions are 'ultra-compact', that 
is, their densities are very high, so that the radius of the 
cavity is very small. Indeed, the estimated gas densities of 
their sample exceed about 10 4 cm -3 (Hunter et al. 2000). 
For such high density region, the Stromgren radius is very 
small, less than 0.1 pc. It is corresponds to < 5" for the 
typical distance of their sample (about 5 kpc). Since the 
cavity's radius is less than the Stromgren radius, the cavity 
cannot be resolved by their observations. If we select some 
nearby spherical 'compact' H n regions, whose gas densi- 
ties are about 10 3 cm~ 3 , as the targets, we may detect the 
dust cavity by SCUBA. 

5. CONCLUSIONS 

We examine a possible distribution of dust grains in H II 
regions. Assuming the geometry with the central dust cav- 
ity, which is suggested theoretically and observationally in 
the literature, we can determine the cavity's radius by us- 
ing a simple transfer model of Lyman continuum photons 
in a dusty spherical medium. We adopt the extinction law 
recently developed by Weingartner & Draine (2001), which 
is based on two components of grains, carbonaceous and 
silicate grains, and is a function of the ratio of visual ex- 
tinction to color excess, Ry. We assume Ry = 5.5, which 
is suitable for high density regions. In the transfer model, 
only one free parameter is included; the radius of the dust 
cavity. 

The fraction of Lyman continuum photons contributing 
to hydrogen ionization (i.e. ionizing photons) is adopted as 
a new constraint to determine the cavity's radius. The ion- 
izing photon fractions of 13 'compact' or 'ultra-compact' 
H II regions in the Galaxy are determined from the ratio 
of infrared to radio fluxes via the method developed by us 
recently (Paper I and II). Its mean is 0.55. It is consistent 
with our previous results; a half of LC photons is absorbed 
by dust before they ionize the neutral hydrogen. 

We examine the effect of various quantities on the re- 
lation between the ionizing photon fraction and the dust 
cavity's radius. Setting a fixed radius of the cavity, we find 
that the ionizing photon fraction decreases as the ioniza- 
tion parameter increases. Also we find that the fraction 
decreases as the value of Ry decreases or as electron tem- 
perature increases. 

We determine the dust cavity's radius of individual sam- 
ple region so as to reproduce the fraction of ionizing pho- 
tons estimated from observational data by the transfer 
model. The mean determined radius of the dust cavity 
is 0.39 times Stromgren radius for all of the sample. We 
divide our sample into two subsamples: one is the 'ultra- 
compact' sample, and the other is the 'compact' sample. 
Our classification is based on the Stromgren (or ionized) 
radius of the sample regions, and is consistent with that 
of Habing & Israel (1979) based on the gas densities. A 
typical radius of the dust cavity for the 'compact' sample 
is 0.30 ± 0.12 in units of the Stromgren radius, where the 



error includes both uncertainties of the individual estimate 
and the variation among sample regions. This corresponds 
to about 40% of the ionized radius and about 0.3 pc. Also, 
a typical cavity's radius for the 'ultra-compact' sample is 
0.48 ± 0.26 in units of the the Stromgren radius, which 
corresponds to about 60% of the ionized radius and about 
0.2 pc. 

Since the uncertainties of these values are somewhat 
large and the 'ultra-compact' sample consists of only three 
regions, it is uncertain whether the difference of the cav- 
ity's radius between two subsamples is real and indicates 
an evolutionary sequence of H n regions. In any case, we 
can reject the dust distribution without the central cavity 
with the significance level about 0.001% for both samples. 
Thus, wc conclude that the model nebula filled with dust 
uniformly is invalid in order to explain the ionizing pho- 
ton fraction. The dust cavity model is supported indepen- 
dently of the previous works. 

We discuss the formation mechanism of the dust cavity. 
Although the radiation pressure and stellar wind by the 
central source can produce the central cavity, it is uncer- 
tain which is dominant mechanism. On the other hand, 
the dust sublimation process is not a major mechanism 
to form the cavity. This is because the expected cavity's 
radius by the sublimation process is too small to explain 
that obtained by us. 

We discuss the detectability of the dust cavity by present 
and future infrared-submillimcter facilities. SIRTF is the 
most powerful facility to detect the cavity. Indeed, we ex- 
pect that SIRTF can detect the cavities in almost all H n 
regions of our Galaxy. SCUBA can also detect the cavities 
if we select nearby (typically < 5 kpc) spherical 'compact' 
H II regions, whose gas densities are about 10 3 cm" 3 , as 
the targets. Moreover, Japanese IR satellite, ASTRO-F 
can detect the cavities in H n regions located within only 
1 kpc from us. Since ASTRO-F surveys the whole of the 
sky, we may detect a lot of closest H II regions with the 
central cavity. 

The author would like to thank anonymous referee for 
many useful comments to improve this paper extensively, 
and H. Kamaya and H. Hirashita for their suggestions and 
discussions to carry out this work. The author has made 
extensive use of NASA's Astrophysics Data System Ab- 
stract Service (ADS). 
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Table 1 

Observational properties of Galactic H ii regions 



Object 


R 


D 


S5GHZ 


-flR 




n c 


Class 




(kpc) 


(kpc) 


(Jy) 


(10~ 14 W cm- 2 ) 


(K) 


(cm -3 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Gl. 13-0.11 


0.2 


8.5 


1.6 


2.1 


5900 


1000 


C 


W31B4 


3.2 


5.5 


2.9 


2.4 


6800 


1000 


C 


G23.95+0.15 


4.3 


5.2 


1.78 


2.3 


6000 


3500 


u 


G29.96-0.02 


4.4 


8.5 


3.74 


8.2 


6100 


1500 


c 


G25.38-0.18 


5.3 


11.5 


2.5 


3.9 


6000 


1600 


c 


G333.60-0.21 


6.0 


3.0 


40 


51 


6200 


4400 


u 


G45.45+0.06 


6.5 


8.3 


2.5 


3.8 


7200 


1000 


c 


W51e 


6.7 


7.3 


11.5 


13 


7100 


1600 


c 


NGC3576 a 


7.9 


3.6 


36.8 


22.5 


7500 


6500 


u 


NGC3603P1 


8.9 


7.2 


6.76 


6.1 


6900 


1100 


c 


NGC3603P2 


8.9 


7.2 


4.14 


4.5 


6900 


1100 


c 


G298. 22-0.34 


9.9 


10.5 


22.2 


7.4 


8600 


2200 




S156 


10.2 


3.4 


1.4 


1.1 


6400 


1000 





a Cols.(4), (5), and (7) are averaged by two different observations. 



Note. — Col. (1): Object name. Col. (2): Galactocentric radius. Col. (3): Distance from 
the Sun. Col. (4): Flux density at 5 GHz in each beam size which is selected to match the 
IR beam. Col. (5): Integrated IR flux in each beam size. The wavelength range is 18-160 
/im. Col. (6): Electron temperature. Col. (7): Electron density from [OIII] doublets (52, 
88 ^m). Col. (8): Classification in section 3. 'C denotes 'compact' H 11 regions, and 'IP 
denotes 'ultra-compact' H 11 regions. All quantities are taken from Simpson et al. (1995). 
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Table 2 



Determined properties of Galactic H ii regions 



Object 




-^IR 


/ 


log 


u 


Rs 


n 


yd 




(io 4 ^- 1 ) 


(10 6 L Q ) 








(PC) 


(PC) 




(1) 


(2) 


(3) 


(4) 


(E 


>) 


(6) 


(7) 


(8) 


Gl. 13-0.11 


1.3 


0.59 


0.44 


-1. 


,94 


0.84 


0.64 


0.18 


W31B4 


0.93 


0.28 


0.63 


-2. 


.07 


0.69 


0.59 


0.38 


G23.95+0.15 


0.54 


0.24 


0.45 


-1 


.89 


0.27 


0.21 


0.24 


G29.96-0.02 


3.0 


2.3 


0.28 


-1. 


.70 


1.00 


0.65 


0.15 


G25.38-0.18 


3.7 


2.0 


0.38 


-1. 


.70 


0.92 


0.67 


0.29 


G333.60-0.21 


4.0 


1.8 


0.45 


-1. 


.58 


0.46 


0.35 


0.45 


G45.45+0.06 


1.8 


1.0 


0.36 


-1. 


.91 


1.05 


0.75 


0.14 


W51e 


6.4 


2.7 


0.47 


-1. 


.70 


1.07 


0.83 


0.44 


NGC3576 


4.8 


1.1 


0.77 


-1. 


.62 


0.33 


0.30 


0.75 


NGC3603P1 


3.7 


1.2 


0.58 


-1. 


,85 


1.06 


0.89 


0.47 


NGC3603P2 


2.3 


0.91 


0.49 


-1. 


.90 


0.95 


0.75 


0.34 


G298. 22-0.34 


23. 


3.2 


1.17 


-1. 


.64 


1.04 


1.10 


1.00 


S156 


0.18 


0.050 


0.67 


-2 


,31 


0.38 


0.34 


0.22 



Note. — Col.(l): Object name. Col. (2): Production rate of ionizing photons 
estimated from 5GHz flux density. Col. (3): Total IR luminosity estimated from 
observed IR flux. The wavelength range is 8-1000 fim. Col. (4): Determined 
fraction of ionizing photons. Mean uncertainty for individual value is estimated 
to be ±0.2. Col. (5): Ionization parameter corrected from Lyman continuum 
extinction. Col. (6): Stromgren radius corrected from Lyman continuum ex- 
tinction in parsec unit. Col. (7): Ionized radius in parsec unit. Col. (8): Dust 
cavity's radius normalized by Rs- The adopted Ry is 5.5. Mean uncertainty 
for individual value is estimated to be ±0.2. 
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Table 3 



Mean properties of 'compact' and 'ultra-compact' H ii regions 



Sample 


Number 


T e 


n c 


N LC 


Rs 


n 


/ 


JJd 


rd/n 








(K) 


(cm- 3 ) 


(10 49 cm" 1 ) 


(PC) 


(PC) 








(PC) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


compact 


8 


6600 


1200 


6.8 


0.95 


0.72 


0.45 


0.30 


0.38 


0.28 






±630 


±400 


±3.9 


±0.13 


±0.098 


±0.10 


±0.12 


±0.14 


±0.13 


ultra-compact 


3 


6600 


4800 


5.4 


0.35 


0.29 


0.56 


0.48 


0.57 


0.17 






±1100 


±1200 


±3.2 


±0.073 


±0.062 


±0.15 


±0.21 


±0.20 


±0.080 



Note. — Col. (2): Number of sample regions. Col. (5): Intrinsic production rate of Lyman continuum photons Col. (7): 
Ionized radius in parsec unit. Col. (10): Ratio of the cavity's radius to ionized radius. Col. (11): Cavity's radius in parsec 
unit. We also show the standard deviations of all quantities. 
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Mean dust cavity's radius for various Rv 



Sample 




Rv 






5.5 


4.0 


3.1 


all 


0.39 


0.50 


0.58 


compact 


0.30 


0.43 


0.52 


ultra-compact 


0.48 


0.57 


0.63 
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Yd 

Fig. 1. — Model calculation of ionizing photon fraction, /, as a function of the dust cavity's radius in unit of the Stromgren radius, y^. 
For all panels, the dotted, solid, and dashed lines represent the cases of the logarithmic ionization parameter, U, of -2.5, -2.0, and -1.5, 
respectively. The electron temperature, T e , is set to be 6600 K. In panels (a), (b), and (c), we show the cases of the ratio of visual extinction 
to color excess, Rv = 3.1, 4.0, and 5.5, respectively. 
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Fig. 2. — Same as Figure 1. In panel (a), we show the variation of the f-yd relation of the case of log U = —2.0 and T e = 6600 K for various 
Rv values. In panel (b), we show the variation of the case of logU = —2.0 and Ry = 5.5 for various T e . 
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Fig. 3. — Histograms of the determined quantities, (a) Ionizing photon fraction, (b) logarithmic ionization parameter, (c) Stromgran radius 
in parscc unit, (d) Ionized radius in parsec unit, (e) Dust cavity's radius in unit of the Stromgrcn radius. The white histograms represents 
the 'compact' H II regions sample, and the shaded histograms represents the 'ultra-compact' sample. The dotted hisograms are two objects 
removed from further discussions. 



